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' We give an elementary introduction to a recent diagrammatic extension of dynamical 

1— ( ' mean field theory (DMFT) coined dynamical vertex approximation (D^TA). This approach 

^ , contains the important local correlations of DMFT, giving, among others, rise to quasipar- 

;_( ■ tide renormalizations, Mott-Hubbard transitions and magnetism, but also non-local corre- 

"^ ' lations beyond. The latter are at the very essence of many physical phenomena in strongly 

; , correlated elecectron systems. As correlations are treated equally on all length scales, DPA 

jrt ■ allows us to describe physical phenomena such as magnons, quantum criticality, and the 

a interplay between antiferromagnetism and superconductivity. We review results hitherto 
obtained for the Hubbard model in dimensions d = 3, 2, and 1. 

rH ■ 

Q ' §1. Introduction 

\ Dynamical mean field theory (DMFT)^' ■*' has been a success story for the the- 

T— I ■ oretical modeling of strongly correlated electrons systems since it contains the major 

'^ , contribution of electronic correlations: the local ones. Diagrammatically, DMFT 

\^ \ corresponds to all topologically distinct Feynman diagrams, albeit only for their lo- 

00 ■ cal part. Hence, it is a non-perturbative approach which, from another point of view, 

^] , can be envisaged as a mean-field theory where the local Coulomb interaction U on 

C^^ I all sites but one is replaced by a dynamical (frequency-dependent) mean field given 

r^ ■ in terms of the self energy E. First, DMFT was applied to model Hamiltonians for 

f^ , studying, among others, the Mott-Hubbard transition,^'' ■^-''^' magnetism, ^)"^'') and 

^ I the strong quasiparticle renormalization of correlated electrons, including the only 

r'"! ■ recently observed kink in the quasiparticle's dispersion relation. ^^'^^-^ In the last 

r^ \ years, the focus shifted towards the realistic calculation of material properties'^' ^^' 

C^ ' and the inclusion of non-local correlations, which are the topic of our paper. 

While usually smaller in absolute terms than the local correlations, the non-local 
correlations nonetheless become essential at low temperatures where they lead to 
pertinent correlation effects such as spin fluctuations, magnons, quantum criticality, 
and (possibly) superconductivity. Actually, many of the most fascinating phenomena 
of strongly correlated electron systems are driven by non-local correlations. 

Two different paths have been followed to extend DMFT by non-local correla- 
tions, see Fig. [TJ Cluster extensions consider instead of the single DMFT site (the 
local problem or Anderson impurity model) a couple of sites and take correlations 
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Fig. 1. Two ways to include non-local correlations beyond DMFT. In cluster extensions of DMFT 
(left), a couple of sites (here two) are embedded in a dynamical mean field which can be formu- 
lated through the self energy S. Restricted to a single site, DMFT is recovered. In diagrammatic 
extensions, the local contribution of DMFT is supplemented by certain non-local contributions. 
Shown is the dynamical vertex approximation (D^TA) where the fully irreducible vertex is purely 
local. These local building blocks are connected by non-local Green function lines dj, result- 
ing in a non-local reducible vertex and self energy. Reducing G,;-, to its local contribution Ga 
reproduces DMFT. 



between these sites into account, see Ref. 19) for a review. There are two main vari- 
ants: cellular (or cluster) DMFT^'')"'^^) where the self energy is directly related to 
a supercell in real space, and the dynamical cluster approximation (DCA)^^''^^' ^^^ 
where it is patched in k space. A generalization in form of a self energy functional 
approach has been formulated more recently by Potthoff.^^''^^^ Also the 1/d ex- 
pansion of OMFT^^''^*^' can be subsumed here since it connects a single-site cluster 
and a two-site cluster. These approaches have been applied successfully for studying 
the one-dimensionaP^^'^^-* and two-dimensional"^^" ^^''^^'^^-^ Hubbard model, as well 
as for the investigation of spin-Peierls physics in Ti203,^^' NaV205,^^^ and V02.^'^^ 
The biggest drawback of the cluster extensions is their restriction to short-range cor- 
relations because the numerical effort is exponentially growing with the number of 
cluster sites, similar to the direct quantum Monte Carlo simulation^^^ of the lattice 
problem. For the one-band Hubbard model, approximately 30-60 sites are possible 
which means a correlation length of about 3 sites in two dimensions and 1 site in 
three dimensions. In realistic multi-orbital calculations, the effort increases further 
with the number of orbitals. Hence, hitherto only two sites have been considered. 

Obviously, we need another route to tackle long-range correlations. This is pos- 
sible through diagrammatic extensions of DMFT which were addressed recently by 
serveral groups. Kuchinskii et al.^^''"^^' combined the local DMFT self energy of the 
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Dynamical vertex approximation {DTK): 

n-particle fully irreducible vertex Fi^ assumed to be local 



n = 1 ^ 1-particle fully irreducible vertex (self energy); yields DMFT 

n = 2 ^ D/^A as discussed in this paper: 

from 2-particIe irreducible vertex Fir 
construct S {local and non-local diagrams) 

n = oo — > exact solution 



Fig. 2. D_rA as the natural extension of DMFT to the n-particle fully irreducible vertex. 



Hubbard model with the non-local self energy of the spin-Fermion model and in- 
cluded long-range correlations this way. Rubtsov et al.^"^' proposed a dual Fermion 
approach to this end. Tokar and Monnier^^' developed a perturbative extension 
which however only includes short-range correlations. We think that a natural ex- 
tension of DMFT coined dynamical vertex approximation {DFK),'^'^' which also in- 
cludes the spin-Fermion diagrams in a systematic way,^^^ is particularly promising. 
Instead of assuming the one-particle fully irreducible vertex, i.e., the self energy, to 
be purely local as in DMFT, Y)FK assumes the same for the n-particle fully irre- 
ducible vertex Fir^'^^' see Fig. [2j These local building blocks are then connected by 
non-local Green functions, yielding a non-local reducible vertex and self-energy; see 
diagrams in the right panel of Fig. [TJ Most natural (and feasible) is the case n = 2, 
i.e., the two-particle vertex. But, in principle, implementations with higher n are 
possible and the exact solution is recovered in the limit n — > oo. The T)FK has been 
proposed by Toschi et al.'^'^' and independently (with minor differences) by Slezak et 
al.,^^' who employ a DC A cluster as a starting point and a perturbative calculation 
of the one-frequency vertex. Very similar ideas including less diagrams have also 
been put forward independently by Kusonose.^^^ 

The paper is organized as follows: Section [2] focuses on the DFA approach and 
algorithm, introducing the full version of the approach in Section 12.11 a restriction 
to the most important diagrams for magnetic fluctuations in Section 12.21 and a 
Moriyaesque A correction in Section 12.31 Physical results for the Hubbard model in 
d = 3, 2, and 1 are presented in Sections 13. H 13.21 and 13. 3^ respectively. Finally, 
Section U] gives a summary and outlook. 
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Fig. 3. Full (self-consistent) UFA algorithm. 



§2. Br A algorithm 

2.1. Full (self-consistent) DFA 

Taking the locality of the fully irreducible*^! vertex as a starting point, we can for- 
mulate the DI^A approach as following, see Fig. [3] for the algorithm: Starting with 
an arbitrary local Green function, we first have to calculate the fully irreducible 
vertex. This is well defined diagrammatically, but in practice one will solve the 
Anderson impurity model numerically using, for example, the exact diagonalization 
(ED)^"^-* or quantum Monte Carlo simulation (QMC) in its Hirsch-Fye,^^^ projec- 
tive, ^^^ continuous time^'^^'^^^ or hybrid^^^ variant. In these numerical approaches, 
one actually does not compute the irreducible vertex directly, but instead the local 
three-frequency (generalized) spin and charge susceptibility and from this -through 
the local parquet equations^^' ^^^- the fully irreducible vertex; see Ref. 44) for details. 

In a second step, the (non-local) reducible vertex is calculated from this (lo- 
cal) fully irreducible vertex. Fig. [1] shows as an example one reducible diagram 
constructed from the irreducible vertex. But of course, we have to sum over all 
reducible Feynman diagrams 



=*) 



This can be achieved more elegantly through the 
parquet equations, ^'^^"^^•* which is still a challenge for an actual implementation of 
this full D/^A approach. 

Having the non-local reducible vertex, one can easily obtain the non-local DI^A 
self energy (third step in Fig. [3|) through an exact equation which follows from the 



*' That is, cutting two Green function lines of the Feynman diagram for the vertex does not 
divide the vertex into two parts. 

**' That is, all Feynman diagrams, which have the local irreducble vertex and the non-local 
Green function as building blocks, and which do not separate into two parts when cutting two 
Green function lines. 
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Fig. 4. Diagrammatic representation of the equation of motion connecting reducible vertex and 
self-energy. 




k"+q,cj 



LXXX> 



;?^x?^ 



tXXA^ 




kt k+q,t 



kT kt k+q,t k+q't kt 



Fig. 5. The two particle-hole channels contributing to the self-energy in the ladder approximation, 
in the longitudinal and transversal spin channel, respectively. Instead of the bare interaction, 
ladder diagrams are constructed from the local vertices (crosshatched) irreducible in the corre- 
sponding spin and charge channels, [reproduced from Ref. 44)] 



Heisenberg equation of motion:^^"^^) 



-q,!^'+cij 



^k' I/' G'k+q,i/H 



(2.1) 



k'q 



see Fig. [J] for the diagrammatic representation. Here, k, k' , and q denote the three 
involved wave vectors; v, u' , and lo the corresponding Matsubara frequency; n the 
number of electrons per site, U the interaction and T the temperature; -T^k'q ^^ ^^^ 
afore calculated reducible four-point vertex and Gk,jy the Fourier-transform of the 
non-local Green function Gij taken as a starting point in the first line of Figures [3l 
From the self energy Eq. (j2-ip . the dispersion relation ek of the lattice and the 
chemical potential ;U, we can now calculate a new local and non-local Green function 
via the Dyson equation 



G\,^u = [ii^n - ek + /" - ^k,uY 



(2-2) 



After this fourth step of the algorithm, we can go back to step 1 and iterate the 
DFA algorithm Fig. [3] until convergence is reached. 

2.2. Restriction to the two-particle-hole channels 

Instead of this full (self-consistent) DFA approach, we have, in a first step, 
implemented a simplified scheme:*^^-* Neglecting the particle-particle channel of the 
Parquet diagrams, we restricted ourselves to the two particle-hole channels which 
are longitudinal and transversal, respectively, in the spin indices; see Fig. [5] b) 
and c). This restriction is reasonable for spin fluctuations where these channels are 
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Fig. 6. Restriction of the D/"A algorithm to particle-hole ladder diagrams. 



dominating and actually diverging close to the antiferromagnetic phase transition 
in two dimensions. For other physical phenomena such as superconductivity the 
particle-particle channel is, as a matter of course, essential. 

The algorithm is shown in Fig. [6l We first calculate the irreducible vertex -Tgrc), ir 
in the spin and charge channel from the numerically obtained susceptibility of the 
Anderson impurity model (step one in Fig. [6l), see Eqs. (9)-(ll) of Ref. 44). From 



this we calculate the reducible vertex FK'^.'^ 

s(c),q 

summation which yields 



(step two in Fig. [6|) through the ladder 
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(2-3) 



where Xoaui — ~^ Sk G'k,i/'G'k+q,jy'+a; is the bubble contribution to the susceptibility. 
In the third step of Fig. [6l the self energy is obtained from these vertices through 
Eq. ([2T]) which yields 
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for the restriction to the ladder diagrams in the two particle-hole channels of Fig. [5] 
b) and c), respectively. Here, -r5(c),ioc is the local (q-summed) counterpart of i~'s(c)^q. 
With the restriction to the two particle-hole channels, actually a non-self consistent 
calculation is preferable so that the algorithm ends with a calculation of the DFA 
Green function from the self energy of Eq. ()2-4p , using the Dyson Eq. (|2-2p . 



2.3. Moriyaesque A correction 

As a numerically inexpensive alternative to the fully self-consistent DFA al- 
gorithm, we recently introduced a Moriyaeque A correction^^^ to the algorithm of 
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Fig. 7. D_rA algorithm with Moriyaesque A correction. 



Section I2.2i In the Moriya theory of weak itinerant magnets, ^^-^ the perturbatively 
obtained susceptibiUty is corrected as follows 



X%^ — ' [{X%^) V A] , 



(2-5) 



where xku — Yli/u' Xsq is the susceptibility. We do the same but for strongly cor- 
related electron systems in the A-corrected D/^A. In our approach, the value of A is 
fixed by the sumrule 
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C/^n(l - n/2)/2. 



(2-6) 



Eq. ()2-5p translates to a A correction of the vertex via 



UU OJ 1—1 



X<i(c),q — [UOqojJ "uu ~ -* s(c),irJ 



(2-7) 



see Ref. 46) for details. Hence in the DFA algorithm, the A correction is one addi- 
tional step, the third step of Fig. [7l Otherwise the algorithm remains unchanged. 

Through the A correction, we include important effects of the self-consistent full 
D/^A. In particular, the reduction of the magnetic transition temperature. In two 
dimensions, even the Mermin- Wagner theorem is fulfilled. ^^^ 
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§3. Results for the Hubbard model 

3.1. Three-dimensional Hubbard model 

Let us turn to the results, starting with the three-dimensional Hubbard model 
on a cubic lattice 

H =-t^cl^Cja + U^ni^n,i. (3-1) 

(ij)cT i 

Here, t denotes the hopping amplitude between nearest-neighbors, U the Coulomb in- 
teraction, cl^{cifj) creates (annihilates) an electron with spin a on site i, riir^ = c-^Cifj . 
In the following, we restrict ourselves to the paramagnetic phase with n = 1 electron 
per site at a finite temperature T. One of the cornerstones of strongly correlated 
electron systems for this paramagnetic phase is the Mott-Hubbard transition, ^^)'^^) 
which was analyzed in DMFT in great detail^-*'^^'^^ with some -by now settled- con- 
troversy for the multi-orbital case.^^^"'^^-' What is less clear however is how far anti- 
ferromagnetic fluctuations above the ordering temperature affect the Mott-Hubbard 
transition. For an orientation, see the phase diagram of the Hubbard model in the 
left inset of Fig. [51 In DMFT, there is no feedback whatsoever of the nearby an- 
tiferromagnetic phase onto the spectrum and the self energy in the paramagnetic 
phase. If we, for example, completely suppress the antiferromagnetic phase through 
"frustrating" the lattice with longer-range hopping, the DMFT spectrum is the very 
same*-*!. D/^A heals this shortcoming of DMFT. 

In our calculations,^'^) we considered three points in the paramagnetic phase of 
the phase diagram, marked as symbols in left inset of Fig. [HI These points are on 
the metallic side of the Mott-Hubbard transition and in the crossover regime where 
the metal continuously transforms into an insulator. In DMFT, the metallic side 
corresponds to a well deflned quasiparticle peak (first point considered at U = ID, 
(3 = 1/T = 0.067Z?, where D = 2v6t is the effective bandwidth/standard deviation 
of the density of states). When increasing U to the crossover regime, this quasi- 
particle peak is first strongly damped, i.e, smeared out (second point considered at 
U = 1.5D, P = 0.0891); for the DMFT spectrum see Fig. [9]). Note, this change from 
a metal to a strongly damped quasi-metal was confirmed recentely by photoemssion 
experiments for V2O3 by Baldassarre et al. ' Upon a further increase of U the 
spectrum develops a pseudogap (third point considered at C/ = 2D, (5 = O.ID). Or, 
if we take vice versa the insulating gap as a starting point, this insulating gap is 
filled, see Ref. 77) for an experimental validation by hands of V2O3. Increasing U 
further, eventually a truly insulating gap is obtained. 

In its main panel. Fig. [8] shows the four-point vertex irreducible in the spin 
channel, which was calculated numerically using exact diagonalization as an impurity 
solver for the Anderson impurity model. ^^^ As we see, the vertex strongly depends on 
all three frequencies (main panel and right inset). It becomes unexpectedly strong if 
an increasing U enhances the electronic correlations. Strong local correlations result 
in a large vetrex which, in the vicinity of an antiferromagnetic phase transition, 



*' except for changes possibly induced by a change of the density of states due the additional 
hopping 
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Fig. 8. Dependenceof the local irreducible vertex _r^ ;^^ ^ on the incoming Fermionic frequency 
v' , for the three different values of U and T indicated as symbols in the left inset, which shows the 
DMFT phase diagram with paramagnetic metallic (PM) , insulating (PI) , and antiferromagnetic 
(AF) phase as well as the crossover region between PM and PI. Right inset: Same as main 
panel but at fixed U — ID and at different uj's (squares: u) = 12nT; triangles: u) = GnT; circles 
u! — 0). [reproduced from Ref. 44)] 



entail strong non-local correlations. 

As discussed in Section 12.21 we can calculate the DFA self energy and fc-resolved 
spectral function from the vertex of Fig. [8l The results for a wave vector k = 
(7r/2,7r/2,7r/2) on the Fermi surface are presented in Fig. [H On the metallic side 
at U = ID, the non-local antiferromagnetic fluctuations lead to a relatively weak 
damping of the quasiparticle peak (QP). Here, the length scales, from which the 
main contributions of the antiferromagnetic fluctuations stem, are relatively long 
range. Physically, we can understand the damping simply from the fact that the 
quasiparticles live in a bath of bosonic spin fluctuations. The quasiparticle can 
scatter at these fluctuations, resulting in a reduced life time, i.e., a larger imaginary 
part of the self energy. 

Entering the crossover regime from the metallic side, the local correlations are 
strongly enhanced and entail strong non-local correlations in the vicinity of the 
antiferromagnetic phase transition. This can be seen for U = 1.5D in Fig. [9l where 
we find a strongly damped quasiparticle peak and a splitting of the Hubbard bands, 
which can be considered to be a precursor of the spin polaron side peaks which 
develop in the antiferromagnetic phase. "^^^ 

Finally, in the more insulating-like part of the crossover regime, aX U = 2D 
the antiferromagnetic correlations become even stronger, albeit they are now more 
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Fig. 9. DFA self energy (left: real part, middle: imaginary part) and spectral function (right) at 
a point k = {n/2,n/2,n/2) on the Fermi surface for U = ID (top), 1.5D (central), and 2D 
(bottom) [blue/light grey line] . In the boxes on the right hand side we note the most important 
differences in comparison with DMFT [red/dark grey line], [reproduced from Ref. 44)] 



short range in nature. The paramagnetic phase becomes hence considerably more 
insulating, indicating a shift of the IVIott-Hubbard transition towards lower values of 
U. It also shows that the IVlott-Hubbard transition is strongly affected by antiferro- 
magnetic fluctuations. 

The results shown here are close to the antiferromagnetic phase transition, where 
the antiferromagnetic fluctuations are particularly strong. If we increase tempera- 
ture, i.e., if we are deep inside the paramagnetic phase, the difference between DFA 
and DIVIFT becomes very small. Already for T > 2Tn, i.e., twice the Neel tempera- 
ture, the differences are minor. 

All results shown were obtained without A correction; and in three dimensions 
the A correction has only a minor (quantitative) effect. It reduces the antiferromag- 
netic transition temperature Tjy since T/v is somewhat overestimated in DIVIFT. But 
the same quantitative effects as in Fig. [9] occur also with A correction albeit at this 
lower T/v, see Ref. 46). 

3.2. Two-dimensional Hubbard model 

Let us turn now to the two-dimensional Hubbard model on a square lattice where 
antiferromagnetic fluctuations are expected to be much more pronounced since the 
particle-hole ladder diagrams diverges logarithmically in two dimensions. Our insight 
into these antiferromagnetic fluctuations has been recently improved not only by the 
cluster extensions of DMFT dealing with short-range correlations, but also by the 
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Fig. 10. Same as in Fig.[^but for the two-dimensional Hubbard model and a single value oiU = D 
at /3 = 17. Shown are results for two different k points on the Fermi surface, in the nodal 
(upper panel) and anti- nodal direction (lower panel). The antiferromagnetic fluctuations lead 
to a complete redistribution of the spectrum in the region of the quasiparticle peak; a pseudogap 
opens. 



functional renormalization group method.^''' ^^^ The latter is very complementary 
to the cluster extensions since fRG deals with correlations on all length-scales, but 
as a (renormalized) perturbative approach it does not include the effects of strong 
correlations such as the Mott-Hubbard transition and strong renormalizations of the 
quasiparticle bands. In D/^A, we cover both of these aspects, strong correlations 
and multiple length scales, in a single approach. 

Let us emphasize that in two-dimensions the Moriyaesque A correction of Section 
12.31 is certainly important to obtain accurate results since otherwise the antiferro- 
magnetic fluctuations will be overestimated. Among others, the antiferromagnetic 
transition temperature T/v would be, without A correction, much too high, i.e., at 
its DMFT mean field value ~ instead of zero as required by the Mermin- Wagner 
theorem. Nonetheless, we will present here results without this A correction which 
can be compared with the reliable A-corrected results of Ref. 46). 

Figure [10] shows the same real and imaginary part of the self energy and spec- 
tral function as in the first line of Figure [9] for U = D, but now for two-dimensions 
(where the effective bandwidth is D = 4t). The difference to three dimensions is 
dramatic since in two-dimensions the antiferromagnetic fiuctuations lead to a com- 
plete redistribution of the spectrum. The quasiparticle peak of DMFT is swallowed 
up in a pseudogap. Even at a relatively weak interaction strength U = D, the ef- 
fects of antiferromagnetic fluctuations are much more pronounced than for a much 
stronger interaction strength in three dimensions. Also the k-dependence is dra- 
matic, in contrast to three dimensions, where we only showed results for one k-point 
because of the weak k-dependence. Instead, in two dimensions, the comparison 
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Fig. 11. Spectral function for the two-dimensional Hubbard model off half-filling in the nodal (left) 
and anti-nodal (right) direction. Parameters are n « 0.8 electrons/site, a next-nearest neighbor 
hopping i' = 0.3t, 13 = 100 and different values of the interaction U = ID, 2D, and 3-D as 
indicated. 



of nodal [k = (7r/2, vr/2)] and anti-nodal direction [k = (vr, 0)] shows a strongly 
anisotropic pseudogap. If we would increase temperatures away from the antifer- 
romagnetic phase transition, the pseudogap would hence first turn metallic in the 
nodal direction. 

Fig. [11] presents results off half-filling for the asymmetric Hubbard model with 
t' = 0.3t in a notation where t and t' have opposite prefactors. At weak coupling, 
U = ID, we clearly see a quasiparticle peak which is already somewhat damped 
by antiferromagnetic fluctuations; Hubbard side bands are hardly visible. Note, 
however the strong anisotropy of the damping which is much stronger in the anti- 
nodal direction. Increasing U to 2D, the damping becomes stronger, but also more 
isotropic. Spectral weight is transferred to the Hubbard side bands. Upon further 
increasing U to 31?, these tendencies continue. The quasiparticle weight is very small; 
and the corresponding quasiparticle peak is strongly damped in an almost identical 
way for the nodal and anti-nodal direction. As we noted already for the three- 
dimensional Hubbard model, the antiferromagnetic fluctuations for weak coupling 
are very long-ranged explaining the pronounced difference between nodal and anti- 
nodal direction. The strong non-local correlations at larger U are somewhat more 
short-range in nature. Consequently, the anisotropy is not so pronounced. 
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Fig. 12. Spectral function of the two-dimensional Hubbard model for different k-points along the 
two high symmetry lines indicated in the upper panel, i.e., in the nodal and anti- nodal direction 
respectively. 



For comparison with angular resolved photoeniission experiments, we also pro- 
vide a scan through the Brillouin zone along lines of high symmetry (in the nodal 
and anti- nodal direction), see upper panel of Fig. [T2j Note, that t' = so that 
the line (0, 0) -^ (vr, 0) is connected to the (shown) line (vr, 0) — > (vr, vr) by particle- 
hole symmetry. When approaching the Fermi surface, the (low-energy) quasiparticle 
peak moves towards the Fermi energy. However, this movement stops at some point 
inducing a pseudogap. Let us note that with the A-correction the antiferromagnetic 
fluctuations are weaker than without A correction^^' but nonetheless still very strong. 
As a consequence, the pseudogap survives only in the anti-nodal direction. 

Most recently, Li et al.^^' calculated also the D/^A lattice susceptibility. The 
results were obtained without A-correction and are similar to those from the dual 
Fermion approach at higher temperatures (and much better than DMFT). How- 
ever, at lower temperatures there are deviations since without A correction DFA 
overestimates antiferromagnetic fluctuations. 

3.3. One- dimensional Hubbard model 

One of the fascinating effects of electronic correlation is the separation of spin 
and charge in one dimension. Slezak et al.^^' addressed this issue in their one- 
dimensional D/^Aesque calculation. As a starting point the authors considered an 
8-site DCA cluster and supplemented it with additional longer range correlation as 
in the DFA, albeit with a restriction of the vertex to the second order diagrams 
shown in Fig. [131 
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Fig. 13. Second order diagrams of the DFAesque calculation by Slezak et al. for the one- 
dimensional Hubbard model, [reproduced from Ref. 46)]. 



This way, the authors where able to go to a much larger number of lattice sites Nc 
than in DC A, see Fig. HM This ahowed Slezak et al. to conclude that for A'c — > oo, 
spin and charge velocity are different so that an electron-like spin-charge excitation 
separates into its spin and charge components with time. This indicates taht the 
inclusion of non-local correlations results in a non-Fermi-liquid phase, suggestively 
a Luttinger liquid. 



§4. Summary and Outlook 



In summary, we have presented the DFA approach in an algorithmic form, 
including the restriction to the two particle-hole channels and the Moriyaesque A- 
correction. Let us also emphasize here that DFA is a complementary approach 
to the extended DMFT (EDMFT).84)-85) EDMFT deals with the local correlations 
induced by non-local interactions, whereas DI^A deals with the non-local correlations 
stemming from a local interaction. The results for the Hubbard model in d = 3, 2, 
and 1 dimension show effects of antiferromagnetic fluctuations such as the damping 
of the quasiparticle peak and the shift of the Mott-Hubbard transition to lower values 
of U in three dimensions, the opening of a pseudo-gap and strong anisotropy of the 
spectrum in two-dimensions, and the spin-charge separation in one dimension. 

Since short- and long-range range antiferromagnetic fluctuations are treated on 
the same footing with a numerical effort much more manageable than for cluster 
extensions of DMFT, we see excellent prospects for this approach in the future. 

On a model level, a variety of fascinating phenomena can be addressed such as 
magnons, quantum criticality and the interplay of antiferromagnetic fluctuations and 
superconductivity. But also realistic multi-orbital calculations are possible for which 
the cluster extensions are too severely restricted because of the numerical effort. Such 
local density approximation (LDA) -|- DI^A calculations would, among others, allow 
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Fig. 14. Spin and charge velocity as a function of system size Nc calculated by DCA (triangles, 
circles) and DFA with a 8-site cluster as a starting point (squares). The parameters for the 
one- dimensional Hubbard model where U — W = 1 {W: bandwidth), k — 7r/2, n = 0.7 
electrons/site, and inverse temperature /3 — 31. [reproduced from Ref. 46)]. 



US to study the effect of antiferromagnetic fluctuations on the Mott-Hubbard transi- 
tion in V2O3 where the model scenario of Section [3. II is realized in an actual material 
and which was previously studied successfully in LDA+DMFT.®^' ^^^ We hope that 
this review, along the lines of a talk given at the Yukawa Institute in November 2007, 
will stimulate more groups to apply the DFA approach to challenging problems of 
strongly correlated electron systems. 
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